Signum Function Method for Generation of Correlated Dichotomic Chains 
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We analyze the signum-generation method for creating random dichotomic sequences with pre- 
scribed correlation properties. The method is based on a binary mapping of the convolution of 
continuous random numbers with some function originated from the Fourier transform of a binary 
correlator. The goal of our study is to reveal conditions under which one can construct binary 
sequences with a given pair correlator. Our results can be used in the construction of superlattices 
and waveguides with selective transport properties. 
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I. INTRODUCTION 

The study of properties of disordered complex systems 
with spatial and/or temporal correlations is one of the 
hot topics in modern physics. Recently, much attention 
was paid to the related problem of how to construct dis- 
ordered materials with specific transport properties that 
are due to underlying correlations in a disorder. One of 
the important applications of this problem is a creation of 
electron nano-devices, optic fibers, rough surfaces, acous- 
tic and electromagnetic waveguides with selective trans- 
port properties. 

It is known (see, for instance, P, 0, [1, 0]) that many 
of properties of systems with weak disorder are de- 
termined by the binary (pair or two-point) correlation 
function of a corresponding random process. Recently 
it was found 0, 0, that specific long-range corre- 
lations in disordered potentials can lead to anomalous 
transport properties. To date, there exist a number 
of algorithms for generating long-range correlated se- 
quences with prescribed correlations. Among such al- 
gorithms the most wi desp read one is the convolution 
method [i, 0, S i m III m, M MMM- In this 
method random elements in the generated chain can be 
of any value from — oo to oo. However, in many appli- 
cations it is more convenient to construct the sequences 
with finite number of random values. An important ex- 
ample of such system, being occurred in nature, is a se- 
quence of nucleotides in a DNA molecule, consisting of 
four elements only. 

The simplest case of a random sequence of finite ele- 
ments is a stochastic dichotomic (binary) chain of only 
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two elements. In contrast with the case of sequences of 
continuous elements, the problem of a construction of bi- 
nary sequences with given correlation properties turns 
out to be tricky. As was recently shown in Ref . (l6l| . 
there is a serious restriction on the type of pair correla- 
tors in binary sequences, in contrast with the sequences 
with continuous distribution of their elements. Many 
of related results for binary sequences can be found in 
Refs.[i3,IIi,[il. 

A direct way to create dichotomic correlated sequences 
is to apply the signum function to the sequence of con- 
tinuous values, obtained with the convolution method. 
This method is based on the convolution of a white- 
noise sequence with some function that is determined 
by the desired pair correlator. However, as was numeri- 
cally found in [20| , the created binary sequence turns out 
to have the pair correlator different from the expected 
one. To date, it remains unclear how to construct binary 
sequences having the same pair correlators as in the se- 
quences with continuous values of their elements. In spite 
of its quite simple form, the signum function method is 
not rigorously analyzed in the literature. 

The present paper makes an effort to clarify this prob- 
lem. Our aim is to understand under what conditions the 
discussed method allows to construct binary sequences 
with a desired pair correlator. We perform a detailed 
theoretical analysis of the restrictions arising for binary 
sequences with long-range correlations. Main attention 
is paid to the step-wise power spectrum resulting in a 
power decay of correlations. This type of correlations is 
extremely important in various applications, such as a 
creation of devices with selective transport properties. 

The paper is organized as follows. Section 2 is devoted 
to general properties of binary sequences. In particu- 
lar, we derive some of the conditions restricting the form 
of pair correlators that a binary sequence can have. In 
Section 3 we describe in details the signum-generation 
method and derive basic relations needed for a further 
analysis. In next Section 4 we analyze the case of a bal- 
anced (unbiased) dichotomic sequence, i.e., the sequence 
with the zero mean-value. Here we also display the re- 
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strictions related to the discussed method. Our main 
findings are reported in Section 5 where we consider the 
most interesting case of long-range correlations with the 
step- wise spectrum. In Section 6 we ask a question about 
general restrictions appearing in the case of power de- 
cay of correlations for binary sequences. In last Section 
we give some additional remarks concerning binary se- 
quences, and summarize our results. The Appendices 
contain some of details of analytical and numerical cal- 
culations. 



II. NECESSARY CONDITIONS 

Let us start with generic properties of dichotomic se- 
quences, not associated with specific choice of a genera- 
tion method. In view of constructing the sequences with 
a given pair correlator, it is of great importance to know 
what are restrictions on the type of correlators allowed 
for dichotomic sequences. To shed light on this problem, 
here we consider a statistically homogeneous random se- 
quence of symbols s„ consisting of the values " — 1" and 



Sn = {-i,i}, nez 



-2,-1,0,1,2,... (1) 



The canonical definition of the correlation function reads 
as follows 

Cs{r) = s;^s;^^s^ ^CsiO)K,ir), (2) 

where s = and Cs (0) = — s'^ are the mean value 
and variance of Sn, respectively. Taking into account a 
peculiar property 



4 = = 1 (3) 
of our dichotomic sequence s„ , one obtains 

C7,(0) = l-s2. (4) 

It should be emphasized that the direct relation between 
the mean value and higher one-point moments is a spe- 
cific property of any dichotomic chain. For instance, for 
the binary chain £„ consisting of "0" and "1" we have the 
relation Ce(0) = £(1 — e) since = £„ in this case. On 
the contrary, for a sequence of continuous random num- 
bers of the Gaussian type the variance and mean value 
are independent parameters . 

To proceed with the two-point moments, we associate 
the correlator Cs (r) with the probabilities of two symbols 
with the same or opposite signs, occurring at the distance 
r, 

4P(±1,^,±1) = (s±l)2 + C,(r), (5a) 

r-l 

4P(±l,^,Tl) = l-s2-C,(r). (5b) 

r-l 



These relationships also are peculiar properties solely of 
dichotomic sequences, they are drastically distinct from 
the analogous relations for other random processes (see, 
e.g., Eq. (|B8[) for the two-point probability density of 
the Gaussian process). This fact is strictly confirmed by 
straightforward calculations of Eqs. ([5|). 

In order to obtain the expression ()5a|) . one should write 



the average (s„ ± l)(s„+r ± 1) via the correlator 



{^n i l)('^n-t-r ^ f ) — ^n^n+r i -^n i ^n+r ~1~ 1 

= Cs{r) + {s±lf. (6) 

On the other hand, the same average can be calculated 
with the use of two-symbol probabilities, 



(s„ + l)(s„+, + 1) = 2 • 2 • P(l,^, 1) 

r-l 

+0-2.P(-l, ... ,l) + 2.0-P(l, ... ,-1) 



-|-0-0-P(-l,_^,-l) = 4P(1,_^, 1). (7) 



(s„-l)(s„+,, -l) = 0-0-P(l,^ 



1) 



+ (-2) • • P(-l,^, 1) + • (-2) • P(l,^, -1) 

r—1 r—1 

+ (-2) • (-2) • P(-l,^, -1) = 4P(-1,^, -1). (8) 

r—1 r—1 

Then, the combination of Eqs. (O, (I?]), and ^ results in 
the equality ([5a| . 

Similarly, in order to derive the expression (|5bp. one 
can write. 



i^n i l)('5n-|-r ~F f) — ^n^n+r "F '^n i ^n-\-r 1 

= Cs{r) + s^-l. (9) 

Again, the above average can be calculated employing 
two-symbol probabilities. Specifically, 



(s„ + l)(s„+, - 1) = 2 • • P(l,^, 1) 

r-l 

+0-0- P(-l,^, 1) + 2 • (-2) • P(l,^, -1) 

r—1 r—1 

+0 • (-2) • P(-l,^, -1) = -4P(1,^, -1). (10) 

r—1 r—1 



(S„ - l)iSn+r + 1) = • 2 • P(l,^, 1) 

r-l 

(-2) • 2 • P(-l,^, 1) + • • P(l,^, -1) 

r—1 r—1 

(-2) • • P(-l,^, -1) = -4P(-1,^, 1). (11) 

r—1 r—1 
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From Eqs. dTO]), and ([TTl) it follows the equality ((5b)) . 

Now, with the use of expressions ([5a|) we express the 
correlation function Cs (r) via the probabilities to occur 
three symbols, 

C4r) + {s±lf = 4 ^ P(±l,^,a, ^ ,±1). 

(12) 

Probability P{ztl,^^^,a, ..^^^ jil) is smaller than 

r' — 1 r—r' — l 

both probabilities P(±l,__^^^^, a) and P(a, ,±1)- 

r' — l 1 — r' — l 

Thus, we can write, 
C,(r) + (s±l)2 

min {P(±l,^, a), P(a, ^ ,±1)} 

= min {1 - s2 _ C,(r'), 1 - - ^(r - /)} 
+ min {(s ± 1)2 + C,(r'), (s ± 1)' + C,(r - r')} 
= 1 - s2 + min { - Cs{r'), -Cs{r - r')} 

+ (s ± 1)2 + min {Cs{r'), C,(r - /)} . (13) 

Here we again have used Eqs. ([5]). Then, according to 
the evident relation 

min{a;, y} + min{— a;, ~y} = — |a; — y\, (14) 

we arrive at the condition 

\Cs{r')~Cs{r-r')\ + Cs{r) ^l-s\ (15) 

Finally, it is convenient to rewrite Eq. (|15p for the nor- 
malized correlator Ks{r), 

\Ks{r')-K,ir-r')\+Ksir) ^1. (16) 

Although we have derived this inequality for < r' < r, a 
simple analysis reveals its validity for arbitrary values of 
r and r' . We would like to stress that the condition p6|) 
is applicable for any mean value s of the dichotomic se- 
quence Sn- However, without a loss of generality in what 
follows we consider binary sequences with the zero mean, 
since the statistical properties of considered sequences do 
not depend on mean values. 

In a similar manner one can obtain second inequality 
with the use of Eq. (|5b| , 

l-s2-a(r)=4 J2 ^(±l,^,a, ^ ,Tl) 

^4 J2 min {P(±l^, a), P(a, ^ ,Tl)} 

= min {1 -s^- Cs{r'), {s T 1)' + C,(r - /)} 

+ min {(s ± 1)2 -I- C,(r'), 1 - - C,{r - r')}. (17) 



Since Cs{r) = Ks(r) in the case of s = 0, Eq. (fT7|) gets a 
simpler form, 

\Ks{r') + Ksir-r')\-Ksir) ^1 for s = 0. (18) 

Thus, applying the necessary conditions (fT6|) and (|17|). 
or (|18|) if s=0, one can identify the functions that can not 
be treated as binary correlators of a dichotomic sequence. 

Inequalities ([16]), pT]) and p8|) have to be met for any 
values of indices r and r' and, therefore, they actually 
represent an infinite set of necessary conditions. Evi- 
dently, in every particular case one should choose the 
strongest condition. On the other hand, Eqs. p^ . (fT7|) 
and (fT5)) are automatically satisfied if one of the indices 
equals zero, r = or r' = 0. The same takes place when 
r = r' . Summarizing all these facts, we can combine 
Eqs. (HH) and (HH) as follows 

max{\Ks{r')±Ks{r-r')\TKsir)} ^ 1, (19) 

r^O, r' 7^ 0, r^r' for s = 0. 

The symbol max{. . .} implies the absolute maximum of 
a function with respect to the indices r, r' . Note that 
Eq. ([TO]) is automatically fulfilled for two limit cases, 
namely, for the delta-correlated (white noise) chain, and 
for the sequences with infinitely long-range correlations. 
As one can see, for both these cases Kg(r) = 5rfi or 
Ks{r) — 1, respectively. 

III. BINARY VERSUS GAUSSIAN 

Now we analyze the construction of a dichotomic se- 
quence 7„ by means of the signum function genera- 
tion (SFG) method. It uses an intermediate correlated 
disorder /3„ obtained as a convolution of the uncorre- 
lated Gaussian noise a„ and modulation function G(n). 
Specifically, the 7-sequence is defined by 

7„ = sign(/3„), (20a) 

OC 

Pn = P+ G(n-n')a„'. (20b) 

n'— — 00 

The initial Gaussian white-noise chain consists of 
stochastic variables a„ with the zero mean, unit variance 
and corresponding probability density. Respectively, 

a = 0, Q!„a„' = Sn,n' , (21a) 

PA{an = a) = ^= exp(-a2/2). (21b) 
V 27r 

The bar over a random symbol or function implies the 
stochastic average. 

Evidently, the constructed dichotomic sequence 7^ 
does not change if one normalizes the intermediate se- 
quence Pn by an arbitrary factor. Hence, without any 
loss of generality, we can admit the variance of /?„ be 
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equal to unity. Then, [3 is the mean value of /?„ and its 
correlation function Kfj{r), 



A>(r) = (/3„-/?)(A 



n+r 



(22) 



is normalized to unity, Kfj{0) = 1. By direct substitution 
of Eq. (|20bp into definition (p2|) . the correlator Knir) is 
readily associated with the modulation function G{n), 



Kpir) 



G{r-n)G{n). 



(23) 



n— — oo 



Owing to evenness of Kplr) — Kp{^r) and in accordance 
with Eq. (|23p. one can also restrict the function G{n) 
to the class of even functions, G{—n) = G{n). Note 
that the condition KpiQ) — 1 gives rise to the following 
normalization for G{n), 



K,{0)= £ G'{n)^l. 



(24) 



n— — oo 



It is convenient to pass to the Fourier transform of 
Eq. (P^l) with the use of the standard expressions. 



1 

2^ 



dk JCi3{k) exp{ikr), (25a) 



/C^(fc)= J2 K0ir)expi-ikr). 



(25b) 



The Fourier transform JCp^k) of the pair correlator Kp{r) 
introduced here, is known as the power spectrum of ran- 
dom /3-chain. Since the correlator Ki3{r) is real and even 
function of r, the power spectrum (j25bp is real, even, 
]Cf3{—k) — JCp{k), and non-negative function of the wave 
number k. Analogously to Eq. (pS)) . for the modulation 
function G{n) one can define its Fourier transform G{k), 
which is also real and even function, G{—k) = Q{k). 
The Fourier representation of Eq. ((23)l reads 



iC0{k)^g^k). 



(26) 



Thus, we arrive at the following expression for the mod- 
ulation function G{n), 



1 



G(n) = — / dk JC^J^ (k) cos{kn) 



(27) 



Evidently, the solution ([27ll automatically satisfies the 
normalization condition (|24p . 

Since the initial chain a„ is a delta-correlated Gaus- 
sian noise, the intermediate variables /3„ also constitute a 
Gaussian random sequence with single-point distribution 
function pb{P) (see Appendix E]) . 



1 



exp 



(28) 



In order to reveal statistical properties of the signum- 
generated dichotomic sequence 7„, one should associate 



its mean value 7, variance C^(0) and pair correlator 
C~f{r) with the corresponding independent characteris- 
tics, namely, the mean value (3 and the modulation func- 
tion G{n) [or, the same, with the intermediate correlator 
Ki3{r)]. According to the definition of average, one can 
express the mean value 7 in terms of (3 via the error 
function plj . 

/>oo 

7 = 7^ = / d/?ps(/3)sign(/3) 



3 /''dxexp(-2;V2) = erf(/3/\/2).(29) 

Similarly to Eq. the variance C^(0) is written in 
the form 



c^(o) ^7^-7' = 1-7' 



(30) 



An important characteristic of the stochastic sequence 
7„ is the correlation function C^{r), 



C-,{r) = t^^t;;^- 72 = G^{Q)K^{r). 



(31) 



Its calculation is performed in Appendix |B1 Here we 
write down only the final equation that relates the cor- 
relator K^{r) to Kp{r), 



(1 - j')K^{r) 



dx 



: exp 



V 1 



(32) 



Note that r.h.s. of the latter equation is not elemen- 
tary function, therefore its analytical study is not simple. 
However, the case (3 = allows one to perform complete 
analytical analysis. 



IV. UNBIASED SEQUENCE 

If the mean value of intermediate chain /3„ vanishes, 
(3 = 0, then due to Eq. ((29)) . the mean value of generated 
sequence 7„ vanishes also, 7 = 0. In this case the relation 
([32]) turns out to be remarkably simplified. 



2 f 

K,{r) = - \ 



2 /•^''('■) dx 



Another equivalent form is 



— arcsin[iir^(r)]. (33) 



Kpir) 



sm 



lK,{r) 



(34) 



From the above relations one can conclude that the 7- 
sequence generated with the discussed signum function 
method, is random. Indeed, the decay of correlations 
with an increase of \r\ in the intermediate /3-chain, also 
results in the decay of correlations in the generated di- 
chotomic 7-sequence. 
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Now it is suitable to rewrite Eq. ([34)1 in the Fourier 
representation, 



IC^{k)=S{K^}{k). 



(35) 



Here the symbol S{-}{k) stands for the operator that 
transforms the function K{r) by the following rule, 

oo 

S{K}{k) = sin \^K{r) exp(-ifcr) (36a) 



OO 

+2^|sin ^K{r) - ^K{r)^ cos{kr). (36b) 



It is important to note that the series (I36a|) can converge 
very slowly. Therefore, in the analytical and numerical 
analysis one has to take into account a lot of terms in 
the sum in order to obtain correct result. To avoid this 
problem, we have used the following trick that is based 
on the second equality (j36bp . Namely, since K{r) — > 
when \r\ —^ oo, the latter sum converges quite rapidly 
according to the asymptotic relation 



7„, if and only if S{K-y}(k) is a non-negative function of 
k. 

In view of the revealed restriction ([M)) . let us make 
a qualitative analysis of Eq. (j36b[) . Note that the first 
summand (1 — 7r/2) is negative. Therefore, if the spec- 
trum IC{k) vanishes within some interval of k and the 
whole function S{K}{k) is non- negative there, conse- 
quently, the third summand is positive, exceeding the 
value (7r/2 — 1). However, in the third summand only 
first few terms determine its sign and give correct esti- 
mate of the magnitude. Thus, one may expect that the 
correlation function of a dichotomic chain, which spec- 
trum is zero in some interval, cannot be generated by 
the SPG method. 

It should be stressed that the condition ([55]) is a quite 
serious restriction to the type of correlators allowed for 
dichotomic sequences. As a demonstration, let us con- 
sider simple example of the one-step additive Markov 
chain of variables e„ = {0, 1} obtained according to the 
conditional probability 



P(e„ = l|e„-i) =e + (e„_i - e) exp(-fcc), (40) 



-K(r) —K^(r), 



|r| ^ cx). (37) 



The substitution of Eq. into Eq. yields the 
following final relation between the modulation func- 
tion G(n) and the correlator (r) of the generating di- 
chotomic noise 7„, 



G(n) = - / dk JS{K. 
Jo ^ 



j}{k) cos(fcn). (38) 



Since the correlator Kj{r) is supposed to be known, the 
relation should be regarded as the expression deter- 
mining the modulation function G{n). 

As one can see, the SFG method for constructing the 
correlated dichotomic sequence 7„ with the zero mean, 
unit variance and prescribed two-point correlator K^{r) 
reduces to the following steps. First, starting from a de- 
sirable profile of K^{r) and employing Eqs. and 
one has to obtain the modulation function G{n). After, 
the correlated sequence 7„ can be generated in accor- 
dance with Eq. ((20l) . However, it is important to take 
into account the restriction that directly follows from 
Eq. (|35p . Specifically, since the power spectrum of any 
random process, in particular /C^(fc), is a non-negative 
function of the wave number k, the function S{K^}{k) 
also has to be non-negative. 



where the parameter kc is the inverse correlation length. 
As is known p^ . this chain has an exponential pair cor- 
relator, ii'cxp(^) — exp ( — fcckl)- It is instructive that 
the function S{Kcyi-p]{k) takes negative values for small 
kc below /c* « 1.099 . . ., see the data in Fig. [TJ Thus, 
the SFG method can reproduce exponential correlator 
only for kc ^ fc*, however, not for small kc in the most 
interesting region of long-range correlations. 



-5C 
I 

X 




-1 





/c=3jt/4 




■ 1 1 1 


k 

c 



S{K^}{k) > for |fc| «C vr. 



(39) 



This condition becomes apparent from Eq. (|38[) , in which 
the function S{Kj}{k) enters as a radicand. In other 
words, with this method the function (r) can be con- 
sidered as a correlator of a dichotomic random sequence 



FIG. 1: (Color online) Dependence S{ exp( — fcc|r|)}(fc) on the 
correlation parameter kc for several values of wave number k. 
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V. LONG-RANGE CORRELATORS WITH 
STEP- WISE SPECTRUM 

A. Maximal Jump 

Here we demonstrate that the discussed method can- 
not be apphed for a construction of dichotomic sequences 
with long-range correlators resulting in the step-wise 
power spectrum 



sin(A:c'') 



(41a) 



-e(fc,-|fc|), Q<k,^Ti, IfcK^lb) 



This kind of correlations is of specific interest in view of 
applications to ID disordered superlattices with a selec- 
tive transport, see, e.g. [1, 0]. Here kc is the correla- 
tion parameter (inverse correlation length) to be speci- 
fied, and 8(a::) implies the Heaviside unit-step function, 
e(x < 0) = and Q[x > 0) = 1. 









k =^ 










5 


*,=^ 
5 












k - — 
5 


Atc 

k. = — 
5 





























1 2 

k 



FIG. 2: (Color online) Dependence S {sm{kcr)/kcr}{k) on k 
for different values of fee. 



within the interval 
= and kc — TT will 



First, we analyze the values of k, 
< kc < TT. Two specific cases of kc 
be considered afterwards. 

In the analysis of the SFG method the crucial charac- 
teristic is the function S{K^}{k) defined by Eq. ([36|) . It 
has to be non-negative for any value of the argument k 
within the interval |A;| ^ tt, see Eq. (p9|. In the case of 
the long-range correlator (|41ap it is suitable to use the 
following explicit expression for S{Ky}{k), 



S 



sm{kcr) 

kcT 

oc 

;5]{sin 



(fc) 



1 



I) 



TT 



Q{kc-\k\) 



TT sin(fccr) 
L2 



kcT 



TT sin(A:cr) 
2 kcr 



I cos(fcr)(42) 



It is remarkable that the summand in the last term be- 
haves as TT^ /24kcr^ when r oo. Hence, at finite kc the 
sum converges quite rapidly and uniformly. Therefore, 
the sum is a continuous function oi k, in particular, at 
k = kc, and can be easily calculated numerically. 

In Fig.[2]the behavior of the radicand (l42|) in Eq. is 
shown for several values of the inverse correlation length 
kc. Since S{-}{k) is an even function of the wave number 
k, the discussion can be restricted by the interval ^ fc ^ 
TT. From Eq. (|^^ and Fig. [2] one can draw the following 
conclusions. 

1 . Due to the last term in expression the function 
S{ sm{kcr) / kcr^ (k) increases with an increase of k 
for all kc within both intervals (0, kc) and (kc, tt). 

2. The negative jump of S { sm{kcr) / kcv} {k) occurs 
at k = kc, at the same point where the power spec- 
trum (|41bp has a jump. The maximal and minimal 
values of the function are achieved at fc = fcc — and 
k = kc + 0, respectively. The jump is exclusively 
formed by the second term in Eq. (|^^ . Therefore, 
its value reads 

S{ sin{kcr) / kcV^ {k — kc — 0) 

-S{ sm{kcr) /kcr] {k ^ kc + 0) ^ TT^/2kc. (43) 

3. The minimal value 5{ sin(A;c'')/^c''}(^ = kc + 0) 
is negative for all finite values of kc within the in- 
terval < kc < TT. The complete dependence of 
positive function —S{sin{kcr)/kcr]{kc + 0) on kc 
is depicted in Fig. [3] 




FIG. 3; (Color online) Dependence — <S{ sm{kcr) / kcr] (kc + O) 
vs kc in log-scale. 

Item |3] displays clearly that for all finite values kc < tt 
the correlator and corresponding power spectrum (|4ip 
cannot be created by making use of the discussed method. 
Indeed, the modulation function G{n) turns out to be of 
complex value, see Eq. (pSjl . 

Now we determine the values of kc for which the func- 
tion (|41ap can not be the pair correlator of a dichotomic 
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sequence regardless of the method of generation. This 
can be done with the use of the results of Sec. [Ill 

First, we demonstrate analytically that this function 
cannot be the correlation function of a dichotomic se- 
quence s„ with arbitrary mean value s for small but fi- 
nite values of kc, namely, for < kc ^ 1. To this end, 
we take Eq. (fT6|) at r' = 1 and r = = [a/kc], where [x] 
is the integer part of x and a is a constant, 



sinfcc sin fee (''a — 1) , sin fcc^^a 

kef a 



Being expanded in 
reads 



kr k^ivg^ 1) k(.Ta 

small parameter fcc, this condition 



( 



cos a 



sma 



)k, + O{kl)^0. (45) 



It is evident that the l.h.s. of Eq. (|45p can be positive 
at some values of a (e.g., for a = 2tt + 7r/4). Thus, the 
requirement (|16p is violated. 

The inequality (|18p can also result in new restrictions 
for allowed values of kc for the case s = 0. Rewriting it 
at r = 2 and r' = 1, we get 



4sinfcr — sin2fcr ^ 2kr. 



(46) 



Numerical analysis shows that this condition does not 
hold true for all finite kc from the interval Q < kc < k* , 
where k* « 2.139 .... It is met only at kc > k* . However, 
since Eq. psp is just necessary condition, one cannot 
guarantee an existence of the correlator (|41aP even at 

At the critical point kc ^ tt the correlator and spec- 
trum (|4T|) reduce to Kj{r) = Sr.o and IC-f{k) — 1. 
This gives rise to the relations, S{Srfl]{k) = 1 and 
G{n) = 6n,o, hence, /3„ = a„. Consequently, one can 
apply the SFG. However, this specific case of fc^, = tt 
is not interesting because from the Gaussian white-noise 
a„-sequence the method fabricates the dichotomic chain 
7„ again of white-noise type. 

As to the singular point kc — 0, here we have Kj{r) = 
1 for the correlator, and /C^(fc) = 27r5(fc) for the power 
spectrum, therefore, the radicand is 5{l}(fc) = 27r^(fc). 
One can see that the SFG is formally applicable. Be- 
sides, the necessary conditions (|45| and ((46|) are satisfied 
automatically. However, this case is a singular one since 
for any arbitrarily small but finite values of kc it is not 
possible to create a dichotomic sequence with the corre- 
lation properties (|4ip neither by the SFG or by any other 
method. 



B. Partial Jump 

Now we extend our analysis to a more general corre- 
lation function that may have various applications. This 
function also results in a step-wise power spectrum, how- 
ever, with an additional parameter h that determines the 



height of step, 



sm{kcr) 

kcT 



(47a) 



/C^,h(fc) = h + {l~h)—<d{kc-\k\)>Q, (47b) 



Q^h^l, < fee ^ TT, |fc| ^ TT. 



Eq. p7)l coincides with Eq. ((4T|) if the step-parameter 
h = 0. Otherwise, when h = 1 the generated 7-sequence 
turns into a dichotomic white noise independently of kc- 
Also, 7„ becomes delta-correlated at kc — tt for arbitrary 
h. Therefore, aA, kc = tt the conclusions of the previous 
subsection are also valid. 

In what follows, it is convenient to analyze finite values 
of kc < TT. The power spectrum (j47bp is an even function 
of the wave number k and has two symmetric jumps at 
the points k = ±kc. For positive < fc ^ tt the spectrum 
abruptly falls down at k = kc from the maximal value 
^-f.h{k < kc) — h + (I — h)'K/kc to the minimal one, 
K^j,h{k > kc) ~ h. Evidently, this jump can be regarded 
as a mobility edge of disordered ID conductors, if the 
maximal value h + {I — h)^ /kc is much larger than the 
minimal one, h. 



1 



l-h n ^ 

h kr 



(48) 



One can see that for finite h well above zero, this is pos- 
sible only for small kc- What is more tricky, for kc ^ \ 
the mobility edge may emerge even in the case when the 
generated 7-sequence is close to a dichotomic white noise, 
i.e., when 1 — /i <C 1. Therefore, one should have. 



0<A;c<l-ft.<l. 



(49) 



The reason of existence of the mobility edge under the 
conditions (|49p is a significant contribution of the second 
term in the coorelator (|47ap . In spite of the fact that it 
has a quite small amplitude 1 — <C 1, this term provides 
extremely long-range correlations with the characteristic 
scale k-^ > (1 - /i)"^ > 1. 

Now we address the function S{^K^ji]{k) that must 
be non-negative in order to construct the correlated se- 
quence 7„ with the use of the SFG method. In accor- 
dance with the definition (j36p . an appropriate analysis 
can be done with the following explicit expression, 

2 

S{K^^um = l-^{l-h) + {l-h)^Q{kc-\k\) 



00 



r=l 



^ sin(M 
2 kcT 



2 ^ ' krr 



I cos(fcr). 



(50) 



As in the previous case (|4ip . the summand in the last 
term behaves as 7r'^/24fcj?r'^ if r ^ 00. Hence, at finite kc 



8 



the sum converges rapidly and uniformly. It is a contin- 
uous function of /c , in particular, at fc = fc^. 

The numerical calculations of 5{i4r^^/i(r)} (fc) per- 
formed for finite < kc < n and ^ A; ^ tt, are shown in 
Fig. m Together with Eq. ((50)) they provide us with the 
following empirical results. 

1. Due to the last term in Eq. (|50p. the function 
S{Kjji^{k) increases with k at arbitrary values of 
kc and h within both intervals (0, kc) and (fcc,7r). 

2. The negative jump of S{K.y^h}ik) occurs at the 
same point k = kc as for the jump of the prescribed 
power spectrum (j47b[) . The maximal and minimal 
values of the function are reached, respectively, at 
fc = fee — and k = fee + for all values of the 
step-parameter h within ^ /i < 1. The jump is 
exclusively related to the third term in Eq. ([50)1 . 
Therefore, its value reads 



S{K^,h}{kc - 0) - S{K^,h]{kc + 0) = (1 - h) — . (51) 

3. Depending on h, the minimum iSji^-y^/i} (fee + 0) 
can be either negative or positive. Also, the value 
oi S {^K^ ^h} [kc + Q) monotonically increases with an 
increase of h. 
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FIG. 4: (Color online) Function S{K^^h{r)}{k) vs k for few 
values of ft at fee = 1.89 (this kc is very close to the minimum 
point of the curve in Fig. [3] and to the minimum point of 
the upper curve in Fig. [SJ . The function is entirely positive if 
h > hmin ■ When h = hmin , the function goes to zero solely at 
one point k = kc + and is positive otherwise. For h < hmin, 
the function is negative within the whole interval kc < fe < tt. 

The presented numerical analysis can be supplemented 
with the following two points. First, from the treatment 
of the case (|4ip , we know that 

S{K-y^h}{kc + 0) <0 for h = 0. (52) 

Second, from the definition ([50)l one can easily reveal that 

S{K^M}{k)^l for h = l. (53) 



Summarizing our results, one can conclude that there 
exists a threshold value of the step-parameter h that we 
refer to as hmin, that separates the region ^ h < hmin 
in which S{K^^!i^{k) has negative values, from the region 
with non- negative values. 



S{K^,h}{k) ;?0 for hr 



It is clear that the threshold hr, 



(see Fig. S]) 



,^n^h^ 1. (54) 
obeys the equation 



(55) 



and depends on the correlation parameter fee- By sub- 
stitution of Eq. ([5(7| into Eq. ([55]) , one can rewrite it in 
explicit form. 



--(1 

2^ 



{sinf^d 

1 

sm{kcr)\ 

> cos[kcr) 



sm(kcr) 



krV 



krV 



(56) 



The numerical solution hmin{kc) of this equation is dis- 
played in Fig. [5] by the upper curve. 

Eq. can be solved analytically at small k^, result- 
ing in 



I (kc 



1 



/96fc, 



1/3 



OikV^) for kc<^l. 



(57) 

This expression exhibits the limit hmiri — > 1 for kc 0. 

Thus, in accordance with the study performed above, 
the dichotomic sequence 7„ with long-range correlator 
and step-wise power spectrum (|47p can be constructed 
by the SFG method only if its parameters kc and h are 
placed onto or above the upper line in Fig.O Only in this 
area of the (fcc, ft.)-plane the condition ([54]) holds true. 
Unfortunately, practically within this whole area the pa- 
rameters kc and h have values of the order of one, and, 
therefore, one cannot satisfy a quite strong requirement 
(|48)) in order to clearly observe a mobility edge. The only 
exception is a narrow vicinity of the point {kc = 0,h = 1). 
Remarkably, in this vicinity due to specific dependence 
([57)1 of hmin(kc), the conditions (|49p can be satisfied, 



0<fcc<l-ft.<l - hmmikc) 



96/cc 



1/3 



< 1, (58) 



and, consequently, the mobility edge can be achieved. So, 
Eq. ([55]) gives us the only (perhaps, just purely theoreti- 
cal) possibility to arrange a mobility edge in the transport 
through the 7-sequence constructed by the SFG method. 

Now we analyze the consequences of the necessary con- 
ditions formulated in Sec. [ill For the long-range cor- 
relator (|47ap the inequalities (fT9|) lead to the following 
restriction with respect to the step-parameter h, 



'i{i?(r,r')} h 



(59) 
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FIG. 5: (Color online) Space of control parameters kc and 
h. Upper curve is the dependence hmin{kc) while low curve 
depicts ho{kc). Within the lowest area in which ^ h < ho, a, 
dichotomic sequence 7„ with the correlator K-y^hif) does not 
exist. The area with hmin ^ h ^ 1 allows for a creation of 
7n with K.y^h{r) by the discussed method. In the intermedi- 
ate region ho ^ h < hmin the SFG method does not work, 
and an existence of a dichotomic sequence with the step-wise 
spectrum remains an open problem. 



at r 7^ 0, r' 7^ and r ^ r' . Here we have introduced the 
function 



sin krj' ^ sin fcc(r 



sin krV 

(60) 



Since by the definition ft, ^ 0, the requirement ([59|) is 
meaningful only if its l.h.s. is positive. Otherwise, it 
is satisfied automatically. The combination of Eq. (jSQ]) 
with the assumption ^ /i ^ 1 gives rise to the relation, 

/lo ^ /i 1, (61) 

where new characteristic quantity ft-o is introduced, 

/lo = 1 -max"^{l,max{i?(r,/)}}. (62) 

This function /lo(^c) is shown in Fig. [5] by low curve. 
Its piecewise shape is caused by the fact that different 
r and r' contribute to within different intervals of 
fee- Finally, when becomes equal or larger than k* 
2.139. . . [see text after Eq. (|46)) ]. we have /iq = and 
the necessary conditions (|6T|) reduce to the initial ones, 

/i ^ 1. 

Thus, taking into account that the required area of the 
parameters fcc and h is determined by the relation, < 
fee < TT, ^ /i ^ 1 of the [kc, ft)-plane, one can summarize 
the following. A dichotomic sequence with long-range 
correlator and step-wise power spectrum (j47l) does not 



exist within the lowest region ^ ft < h^ik^. When 
hoikc) ^ ft < ftmm(fcc), the dichotomic chain cannot 
be created by the SFG method and there is no answer 
whether it can be created by any other method. Finally, 
within the highest zone with hmin{kc) ^ ft ^ 1, one can 
construct desired dichotomic sequences with the use of 
the discussed method. 



C. Predefined Intermediate Spectrum 

With the SFG method, we first generate intermedi- 
ate Gaussian sequence /?„ and after, the dichotomic se- 
quence 7„. Above, we specified the correlator K^ir) 
of a final dichotomic 7-sequence and analyzed the spec- 
trum K,f3{k) = S{K^}(k) of the intermediate Gaussian 
/3-chain, keeping in mind that the latter must be non- 
negative, see Eqs. ([55]) . ([55]) . 

Below, we ask question about the type of the spectrum 
of 7„, that emerges if the intermediate Gaussian sequence 
[3n is assumed to have given pair correlator with step- wise 
power spectrum of the following form. 



Kp{r) 



sm[kcr) 



(63a) 



|fc|) , < /Cc ^ TT, |fc| < 7<.63b) 



jcp{k) = -e(fc. 



In accordance with the relation ([33]) and Fourier trans- 
forms ([^ , the corresponding correlator and power spec- 
trum of the dichotomic 7-sequence read. 



arcsm 



■ arcsm 



sin(fccr) 



sin(fcc'') 

krV 



(64a) 



exp(— jfcr) 



kc 



eikc-\k\) 



arcsm 



r=l 



sin(fccr)i sin(fccr) 



kr 



kc 



|cos(fcr). (64b) 



Here, last representation for the spectrum )C-y{k) is 
similar to that we have employed for the function 5{-}(fc) 
in its analysis [compare with Eqs. ([5^ and ([H])]. It is 
noteworthy to emphasize that the summand in the last 
term of this representation behaves as 2/37rA;^r^ when 
r — > 00. Hence, as above, the sum converges rapidly and 
uniformly. Therefore, it can be easily calculated numer- 
ically. 

Fig. [6] presents the behavior of the power spectrum 
(j64bp for dichotomic 7,1 in the case of the predefined 
step-wise profile ()63bp for the spectrum of intermediate 
sequence /?„. Since lC^{k) is an even function of the wave 
number fc, the presentation is sufficient within ^ /c ^ tt. 
Assuming kc < tt, from Eq. (j64b|) and Fig.[n]one can con- 
clude: 
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Here we have introduced 




FIG. 6: (Color online) Power spectrum K.~,{k) vs k for differ- 
ent values of the correlation parameter fee- 



1. The function (|64bp is positive within the whole in- 
terval ^ TT. Therefore, it truly serves as a power 
spectrum and its inverse Fourier transform (|64aP is 
valid correlator of the generated sequence 7„; 



2. Due to the second and last terms in Eq. (j64bp the 
spectrum /C^(fc) decreases with an increase of k for 
all kc within whole interval (0,7r); 



The negative jump of /C^(fc) occurs at A: = fcc, at 
the same point as for the jump of prescribed inter- 
mediate power spectrum /C^(fc). The maximal and 
minimal values of /C^(fc) at the jump are defined 
at fc = fee — and k — kc + Q, respectively. The 
jump is exclusively formed by the second term in 
Eq. (I64bp . Therefore, its value reads 



/C^(fcc - 0) - /C^(fcc + 0) = 2/kc 



(65) 



4. The value JC^{k = fee — 0) to the left from the jump 
is always positive. The value IC-y{k = kc + 0) to its 
right is also positive for all kc- Their ratio lC^{kc — 
Q)/K-y{kc + 0) being of the order of one for finite 
fee, seems to slightly increase with a decrease of kc- 
However, it is saturated if the inverse correlation 
length kc vanishes. Indeed, from Eq. (j64bp one can 
easily get, 

JC^jkc - 0) _ TT + I{kc) + (7r/2 - l)fc, 
K^{kc + 0) " 



I{kc) + (7r/2 - l)fc, 
7r + /(0) 

m 



(66a) 

6.0096 ... if fcc -> Oj66b) 



I{kc) = 2k, 



Si 

r=l 

sin(fccr) 



arcsm 



sin(fcc'') 



I cos(fccr); 



m 



2 j dx { arcsin 





sma; 



(67a) 



cos a;(67b) 



In spite of the divergence of the jump the 
convergence of the ratio (|66|) at fcc ^ occurs be- 
cause the values of the spectrum /C^(fc) to the left, 
JC^{kc — 0), and to the right, ICj{kc + 0), from the 
jump, increase with decreasing of kc exactly in the 
same manner, K.^{kc — 0) 2[7r + /(0)]/7rfcc and 

/C^(fcc + 0)R:!2/(0)/7rfcc. 

In addition, it can be analytically shown that in the 
limit fee — > 0, the correlator K'^{r) tends to unity, while 
the spectrum K,^{k) turns into the Dirac delta-function. 



lim K~Ar) — — lim arcsin 
lim /C^(fc) = 2Ti5{k). 



sin(fcc'') 



(68b) 



Therefore, as the correlation parameter kc vanishes, the 
final dichotomic 7-sequence becomes to have extremely 
long-range correlations. 

On the contrary, for kc = tt the correlator Kj{r) re- 
duces to the Kronecker delta-symbol, whereas the power 
spectrum lC^{k) degenerates into unity, 

2 

K^{r) ~ — arcsin [J^.o] — ^r,Q (69a) 

TT 

ICj{k) = 1 for kc = TT. (69b) 

Thus, the final dichotomic 7„ reduces to the white-noise 
chain. 

In summary, in spite of the fact that the jump ratio 
ICj (kc — O) /ICj (fee + 0) in the power spectrum /C^ (k) is of 
the order of unity for all kc < tt, we hope that the mo- 
bility edge can be observed at small enough correlation 
parameter kc 1, due to specific delta-function behavior 
ofJCjik) itself [see Fig. IS] and Eq. ((68bll ]. 



VI. POWER CORRELATORS 

A. Power Correlator for Dichotomic Sequence 

Here we consider an important problem of construct- 
ing a sequence with the power correlation function and 
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corresponding spectrum, 

K^A^) = 5r,o + {kc\r\)-P{l-Sr,o), (70a) 

+ i^Up[exp{ik)] + Lip [ exp(-iA:)] ^TOb) 
p>0, kc^l, \k\f^Tr. 

Here p and fc^ are positive real numbers characterizing 
how fast the correlator decreases. Note that the param- 
eter kc cannot be less than one since Kj^p{r) ^ 1. The 
Fourier transform ICj^p{k) of this correlator is expressed 
via the polylogarithm function hiq(z) that is defined by 




z 



(71) 



For K^,p(r) to be the correlator of a stochastic process, 
it is necessary to have /C-y^p(fc) ^ for all fc. This condi- 
tion is satisfied if and only if the following inequality is 
fulfilled, 



fc, ^ [-2Lip(-l)] 



(72) 



This result is due to the shape of spectrum /C^^p(fc) that 
monotonously decreases with an increase of k within the 
interval (0,7r), and reaches its minimal value at fc = tt. 
The r.h.s. of the condition ([7^ can be calculated in limit 
cases, 



[-2Lip(-i; 



7r/2 — cp , P ^ 1, 
l+p"Mn2, p>l, 
c= 1.029.... 



(73) 



In Fig. [7] the area where power spectrum (j70bp is non- 
negative, is located above the dotted lowest curve. 

If K-y^p{r) is the correlator of a dichotomic ran- 
dom sequence generated by the SFG method, then 
S{Kj^p}{k) ^ for all values of fc. Since this function, as 
well as /C^.p(fc), monotonously decreases, the only condi- 
tion is required. 



S{K^,p}{tt) ^ 0. 



(74) 



Expanding in Eq. (|36ap the sin-function into series, we 
get useful expression 



S{K^,p}{k)^l+J2 



1=0 



22'+i(2; + 1)! 



fc-p(2i + l) 



x{Lip(2;+i) [exp(i/c)] + Lip(2;+i) [exp(-ifc)] |.(75) 
For p oo it can be approximately calculated as 

S{K-yAik) ^ 1 + 2sin {nk^P/2) cos/c. (76) 



FIG. 7: (Color online) Various borders of the parameters 
p and fee above which the following relations are fulfilled: 
(a) Eq. ([72} (dotted curve), (b) Eq. (solid curve), (c) 

Eq. (113) (dashed curve), (d) Eq. ([M} (dash-dotted curve). 



Therefore, taking into account that k^ ^ 1, one obtains 
the following asymptotic of the condition ([7H) 



fcc>l+P ^ln3 for p>l. 



(77) 



The area where the discussed SFG method is applicable 
located in Fig. [7] above solid curve. 

Now we analyze the necessary condition for existence 
of the power correlator (|70ap for dichotomic sequence re- 
gardless of the generation method. For the sake of sim- 
plicity we consider the case 7 = 0. From the inequali- 
ties ([ro]) one can obtain the following relation, 

max{\\r'\-P±\r-r'\-P\ T \r\-P} ^ fcP, (78) 

r ^ 0, 0, r' . 

It is easy to see that the maximum with respect to r', 
occurs at r' = ±1 or r' = r ± 1. Therefore, Eq. ([78]) can 
be rewritten as 



max ( 1 ± (r 

r>0 



1)-Pzpr-P} ^ kP. 



The last condition is equivalent to 



kr^ 2 



>i/p 



2-2pln 2, 
In 2, 



1 



P 



P< 1, 
p > 1. 



(79) 



(80) 



In Fig. [7] the border corresponding to this condition is 
designated by dashed curve. It should be noted that 
Eq. ([50)1 is the necessary condition, thus, if it is met, it 
is still not clear whether the correlator with such values 
of parameters p and kc exists. 



B. Predefined Intermediate Power Correlator 

We have found that the dichotomic sequence 7„ with 
the power correlator (|70ap can be constructed by the SFG 
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method for some values of parameters p and fee that meet 
the condition ([7^ . Now let us take the intermediate 
Gaussian sequence /3„ prescribed to have the power cor- 
relator and corresponding spectrum 

Ki3.,p{r) = Sr,o + ikc\r\)-P{l-5r^o), (81a) 
^t},p{k) = 1 

+ fc;r''|Lip[exp(iA:)] + Lip[exp(-ifc)] |$lb) 

P>0, fee ^ 1, |fc| ^ TT. 



Evidently, the condition ([72]) is implied to be met. In 
accordance with the relation ()33p and Fourier transforms 
([25|l , the correlator and power spectrum of the generated 
dichotomic 7-sequence are described as 



Kj„{r) = Sr.o H — arcsm 

TT 



m[ik,\r\)-P]{l-Sr,o), (82a) 



r=l 



{Kr)~P cos(fcr). (82b) 



We can assert that since Eq. (|72l) is satisfied, i.e. the 
spectrum (jSlbp of the intermediate /3-sequence is non- 
negative, then the spectrum (|82b[) of the generated di- 
chotomic 7„ is also non-negative. 

When \r\ 00, the correlator (|82ap tends to zero in 
accordance with the following asymptotic 



K,,pir) ^ (Klrl)-", K - fcc(V2)'/^ 



(83) 



As was shown, the allowed values of kc are expressed 
by Eq. (|72p. Therefore, the scaling parameter fc^ should 
satisfy the condition 



fc^^ [-^Lip(-l)] 



i/p 



(84) 



This condition for possible values of k'^ and p is met in 
the area above the dash-dotted curve in Fig. [71 

Thus, the mapping of the gaussian sequence with the 
power correlation function (j81ap into the binary sequence 
result in the same power for the decrease of the final cor- 
relator (|82ap expressed by Eq. ([55)1 . However, such a be- 
havior of the final correlator occurs only asymptotically, 
for sufficiently large values of |r|. 



VII. CONCLUSION 

First, we would like to emphasize the following point 
that was briefly mentioned in the beginning. Our study 
of the correlation properties of a random dichotomic se- 
quence 7„ is based on the example ([T|) in which two el- 
ements are " — 1" and "1" . On the other hand, there 
is a simple correspondence between this chain and a di- 
chotomic sequence s(n) consisting of two arbitrary sym- 
bols eo and £1, 



e{n) = {eo,£i}, 



neZ 



-2,-1,0,1,2,. 



The correspondence is expressed by the linear relation- 
ship. 



e{n) 



Eq + El Eq — E\ 

T 7n ■ 



(86) 



The choice of the sign is not important. It derermines 
only into what symbol, £0 or £ii the initial values " — 1" 
and "1" are converted. 

In accordance with Eq. (j86p and due to specific prop- 
erties ([5(7)) . ^ of the 7-sequence, the connection between 
the mean values and variances is as follows. 



\n) = 



_ £0 + £1 £0 ~ El _ 

^= ^ T ;;— 7; 



(87a) 



(87b) 



a(0) = £2(n)-£^ 



(£0 - £1)' 



C^(0). (87c) 



Analogously, the two-point correlation function Cg (r) of 
the £-chain is associated with the binary correlation func- 
tion C^(r) of the sequence 7„ as follows 



C,(r) =£(n)£(n + r)-£2 = l£5L^il!c.,(r) , 



(85) 



The comparison of Eqs. (|87cp and ([55)1 makes obvious 
the equality between the normalized correlators K^{r) 
and K^[r\ 

K,{r) = C,{r)/CM = C^{r)/C^{0) = A%(r). (89) 

Thus, our analysis is valid for any dichotomic sequence. 

Our results can be summarized as follows. We have 
shown that the statistical properties of random di- 
chotomic sequences are principally different from those 
known for sequences with a continuous distribution of 
their elements. We were able to find analytically the 
conditions that can be used to know whether a bi- 
nary sequence can have the desired pair correlator. Note 
that these two conditions are necessary only. 

Another important restriction is due to the inequality 
(|39)) derived under quite general assumptions. We have 
shown that even in the well known case of an exponential 
decay of correlations, there are no binary sequences that 
can be created with the SFG method, unless the decay 
is sufficiently strong. This fact is very important in view 
of many applications. 

Our specific interest was in a possibility to create, with 
the considered method, the binary sequences with long- 
range correlations described by Eqs. (|41aP and (j41bp . We 
have analytically found that the function (|41ap can not 
be a pair correlator of any binary sequence. We have 
also examined a more general correlation function [see 
Eq. (|47)) ] that corresponds to the generalization of the 
step-wise power spectrum. Our extensive examination of 
the signum-function method, applied to this correlation 
function, has revealed the regions of parameters kc and h 
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for which the pair correlator can emerge in a binary se- 
quence. Correspondingly, we identified the regions where 
such a pair correlator can not appear in binary sequences. 
For other values of the control parameters we can not give 
definite answer, therefore, a further study is needed. 

Finally, we analyzed an important case of the power 
decay of the pair correlator. Recently, the problem of the 
generation of random processes with power correlations 
has attracted much attention in the literature. Analyzing 
such correlators, we have found that the SFG method in 
principle allows to construct binary sequences with these 
correlators, however, with some restrictions on the values 
of parameters in Eq. ((70|) . 
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APPENDIX B: PAIR CORRELATOR OF 7„ 

Let us derive the pair correlator ^nln+r- Employing 
the standard integral presentation for the signum func- 
tion. 



sign(z) 



dx 



sin(za;) 



(Bl) 



and Eq. (|A4[) . we arrive, in a manner similar to the calcu- 
lation of the characteristic function (ps (t) in Appendix [Xl 
at the expression, 



2 r°° dxi r°° dx2 



7n7n+r = J{Kp{r),/3) 

\ —) 

sexp [sK(3{r)xiX2] cos [(3{xi — SX2)] (B2) 



9 I I ■ exp I 

Jo Jo 



APPENDIX A: PROBABILITY DENSITY OF /3„ 

The standard way to derive the probability density 
Pb{P) of the random process /3„ is due to its charac- 
teristic function (psit) defined by 



ipB{t) = exp[itf3„ 



d(3 pB{[3)eyiT>{itl3). (Al) 



From the last equality in this definition it immediately 
follows that the probability density psiP) is the Fourier 
transform of ipsit), 



1 

2^ 



dtpB{t)exp{~it(i). (A2) 



To start with, we substitute the explicit expression 
(|20bp for /3„ into the definition (|A1[) for characteristic 
function ipB{t). Then, we rewrite the result as an infinite 
product of exponential functions and take into account 
the statistical independence of uncorrelated random vari- 
ables a„. This procedure yields 



exp[it(3n] = exp(i</3) exp [itG(n — n')an'] ■ (A3) 



n — — 00 



In accordance with the Gaussian distribution (|21bp of 
an, its characteristic function is 



ipA(j) = eKpiiran) = / da /9^(q;) exp(ira) 

J —OQ 

= exp(-rV2). (A4) 

The use of Eqs. jXSj, (|^ with r = tG{n - n'), and 
the normalization condition (I24p results in 

(A5) 



cpBit) ^exp{ipt-t^/2). 



As is known, this characteristic function corresponds to 
the Gaussian probability density One can confirm 

this fact by a direct evaluation of the integral in Eq. (IA2[) . 



Thus, we have reduced the problem to the derivation of 
J{K, f3). To solve it, we obtain the derivative of J^{K, /3) 
with respect to K. After some simplifications one gets, 



d 



1 



3 foo ,2,2, 



X exp [ii'a;ia;2] exp [itp{xi — X2)] ■ 

t=-i,i 

To proceed, we write down the following relation that is 
valid for arbitrary real quantities a and b, 

/oo 
dx exp(— a;''/2) exp(aa;) exp{ibx) 
-00 

= \/2^exp [(a^ - 62)/2] exp(m5). (B3) 



Using Eq. (|B3p with a = Kxi and b = —t[3 we inte- 
grate over X2 and make further simplifications. 



_d_ 



:J{K,P)^^J dxiexp[-xlil- K^)/2] 
X exp(-/3V2) exp [i(3{l — K)xi\ . 
Now we change the integration variable xi. 



x\ = xi\/l - if 2. 
Then, applying Eq. (|B3| with 



a = 0, 



1 - K 
1 + 



(B4) 



(B5) 



we perform the integration over xi that gives rise to the 
expression 



exp 



(-T 



/32 



(B6) 
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The general solution of Eq. (|B6p is 



J{K,(3) = J{Q,P) 



K 



TT Jo \/l - a; 



expf-^). (B7) 



1 + x 



It should be noted that Eq. (|B7p can be also obtained 
by means of the two-point probability density that for 
the correlated Gaussian sequence (in with the correlator 
Kj3{r) is defined by 



PB{Pn=P.f3n+r=P') 

: exp ^ - 



(B8) 



1 f /32+/3'2-2if^(r)/3/3'' 



2^^1- Klir) I 2[1-X2(.)] y 

The last step we should take, is to calculate J {0,(3). 



It can be directly obtained from Eq. (|B2p . 



72.(B9) 



This result can be easily explained. Indeed, the condition 
Ki3(r) = implies that the correlations between /3„ and 
(Sn+r disappear, hence, the correlations between j„ and 
jn+r are absent as well. 

As a result of these calculations, we finally get 



^(A>(r),/3)=7' + 



2 /■^''('■) dx 



: exp 



(-T 
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This expression provides Eq. ((32|) . 
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